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ABSTRACT 


In this study, optimal vapor release rate (or pressure) histories have been generated for an 
industrial semi-batch nylon 6 reactor using Pontryagin’s minimum principle. The batch time 
has been taken as the objective function, which is to be minimized. The pressure is constrained 
to lie between a lower and an upper limit. The temperature, a state variable, is also constrained 
to lie between 220'^ C and 270° C in order to ensure single-phase polymerization. Optimization 
has been carried out with a single end-point constraint (on monomer conversion) and a stopping 
condition (obtaining a product having a desired degree of polymerization, fj-n). Techniques have 
been developed to overcome the discontinuities present in the model, as well as to take care 
of state variable constraints. The effects of various physical and computational variables on 
the optimal pressure history and the corresponding batch time have been studied. It is found 
that the optimal batch time is almost 50% of the industrial value used currently. Interestingly, 
the optimal pressure history is quite similar qualitatively with the current practice though 
quantitatively there is a significant difference. Improvements in reactor operation along these 
lines have been reported. 
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CHAPTER 1 


Introduction 

A significant amount of research on the simulation of nylon 6 polymerization has been 
reported in the past few years. These have formed the subject of several reviews [T5]. The 
effects of various operating variables (initial water concentration, temperature, monofunctional 
acid stabilizers, initial monomer concentration, etc.) on the characteristics of the product [viz., 
monomer conversion, degree of polymerization (DP or /z„), polydispersity index (PDI), etc. ] 
are now well established. The early reports have been followed by several studies [6-16] on 
the optimization of nylon 6 reactors. Hoftyzer et al. [6] used dynamic programming to de- 
termine the optimal temperature and water concentration (related to total applied pressure) 
histories, T(t) and [W](t), required for producing nylon 6 having a desired value, of the 
degree of polymerization in the shortest possible reaction time, tj. Unfortunately, Hoftyzer 
et al. presented only semiquantitative results, because of proprietary reasons. Reimschues- 
sel and Nagasubramanian [7,8] optimized combinations of two-stage isothermal reactors, and 
presented quantitative results. These two studies made use of a kinetic scheme involving only 
the major reactions (ring opening, polycondensation, and polyaddition). Naudin ten Cate [9] 
and Mochizuki and Ito [10] used the minimum principle to find optimal temperature profiles in 
tubular reactors. Naudin ten Cate optimized a two stage reactor, the first one at temperature 
To, having a residence time to, in which water is vaporized continuously. This was followed by 
a nonvaporizing tubular reactor having a residence time tj - to- Mochizuki and Ito, however, 
optimized a single reactor. Both these groups of workers used semiempirical approximations for 
the rates of formation of the undesirable cyclic oligomers, and attempted to minimize the for- 
mation of the water extractable compounds (cyclic oligomers, water and unreacted monomer) 
while simultaneously keeping some product characteristics fixed (end-point constraints). These 
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workers assumed no evaporation of water and/or caprolactam and obtained optimal temper- 
ature histories. Again, only semiquantitative results were presented. Dynamic optimization 
of batch (or tubular) reactors using a variety of objective functions and end-point constraints 
has been carried out by our group [11-16] using Pontryagin’s minimum principle [17,18]. These 
studies have shown that the optimal histories are largely dependent, both qualitatively and 
quantitatively, on the choice of the objective function and the constraints. This observation is 
consistent with similar observations made by Ray and Szekely [17] and Denn [18], who have 
compiled the results of several optimization studies for nonpolymeric systems. Most of these 
studies give valuable insight into the optimal operation of ideal reactors, and indeed, these must 
have been used in the design and improvement of nylon 6 reactors in industry over the years. 
Work on the modelling of actual industrial reactors has only just started to appear in the open 
literature [19]. In this study, the operation of one such reactor (semi-batch nylon 6 reactor) has 
been optimized using Pontryagin’s minimum principle [17,18,20,21]. Optimal pressure histories 
have been obtained for the manufacture of three grades of polymer. Substantial improvements 
are possible in the performance of this reactor once these are implemented. 

The industrial reactor studied here is shown schematically in Fig.l. It is a jacketed vessel 
with a low speed anchor or ribbon agitator used for mixing the highly viscous reaction mass. 
The reaction mass (caprolactam, water and other inert additives like Ti02, etc.)is heated by 
condensing vapors in a jacket. As the polymerization takes place the temperature inside the 
reactor goes above 220°C and vaporization of caprolactam and water takes place, resulting 
in build up of pressure. The pressure inside the reactor, p(t), is maintained to conform to a 
desired history by manipulation of a control valve which allows the vapor mixture of nitrogen, 
monomer and water to pass to a condenser. The pressure histories used currently (refered to 
as ‘reference’ conditions) for three different industrial runs producing different grades of nylon 
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6, are shown in Fig. 2 (curves marked ‘ref’). In this figure, the pressure and the reaction time, 
t, have been nondimensionalized using [19] (see nomenclature) 

^ ~ (P Po)/(Pmax,re/ Po) (1) 

T=tltj^ref ( 2 ) 

The values of the maximum pressure, Pmax.re/, and the total reaction time, used currently, 

are not being reported for proprietary reasons. The jacket temperature. T.. is kept constant 
throughout the reaction at the current value, Tj,re/- 
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CHAPTER 2 


FORMULATION 


The kinetic scheme for nylon 6 polymerization is given in Table 1 [4]. This kinetic scheme 
incorporates the three important reactions viz. ring opening, polycondensation and polyad- 
dition, as well as reactions with the cyclic dimer. Because of the unavailability of precise 
information about the rates of the reactions involving higher cychc oligomers, these have not 
been incorporated in Table 1. Since the cyclic dimer constitutes a major share of the total 
cyclic oligomers in the reaction mass, this assumption is justified. The equations for the fifteen 
state variables, x,(mass and energy balance, and moment equations) are given in appendix 1 
[19]. In general the state variable equations can be expressed as: 

^ = /,(x,u);!; = 1,2,---,15 (3) 

where x and u are the vectors of the state and control variables. 

The ODEs (ordinary differential equations) in Eq.3 can be integrated using the D02EJF 
subroutine of the NAG library. This uses Gear’s technique for integrating sets of stiff ODEs 
[22]. The presence of some discontinuities and the stiffness of ODEs make it difficult to use a 
constant error tolerance (TOL) in the computer code. Provision was made in the algorithm for 
self- adjustment of the error tolerance between 10"“* and 10"^^. When D02EJF failed to solve 
the ODEs with the given error tolerance, it came out of the subroutine giving an error message, 
IFAIL=2. At this point a backward step was taken and TOL reduced by a factor of 10. The 
ODEs were again integrated in the forward direction. This backward movement was continued 
with a reduction in TOL, and the integration repeated till the solution was obtained. A similar 
procedure was followed for the backward integration of the adjoint variables. The change in 
viscosity of the reaction mixture and its effect on the mass transfer rates have been incorporated 
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Table 1 

KINETIC SCHEME FOR NYLON 6 POLYMERIZATION 

l.Ring Opening 


2. Polycondensation 


Ci + W _ 5i 

k[=^ 
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s„ + s 
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Sm+n + 11* ]n,m = 1: 2, ■ - • 



to 

f<2 
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•S'n + C*! 5n+i ; n. — 1, 2, ■ • • 


k'--^ 


4. Ring Opening Of Cyclic Dimer 


^4 

C2 + W 52 
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Rate and Equilibrium Constants 


ki = Atexp(-Et/Rr) + A‘exp{-EURT)Y,([Sn]) 


n=l 



Ki 

= exp [(A 5,- 

-AHifT) /R] 

, i — 1, 2, • • • 5 5 



i 


Ef 



AH, 



(kg/mol-h) 

(cal/mol) 

(kg^/moF-h) 

(cal/mol) 

(cal/mol) 

(cal/mol-K) 

1 

5.9874 X 10^ 

1.9880 X 10'* 

4.3075 X 10^ 

1.8806 X 10^ 

+1.9180 X 10® 

-7.88460 

2 

1.8942 X 10^° 

2.3271 X lO-* 

1.2114 X 10^° 

2.0670 X 10“ 

-5.9458 X 10® 

+0.94374 

3 

2.8558 X 10^ 

2.2845 X 10“* 

1.6377 X 10^° 

2.0107 X 10“ 

-4.0438 X 10® 

-6.94570 

4 

8.5778 X 

4.2000 X 10^ 

2.3307 X 10^^ 

3.7400 X 10“ 

-9.6000 X 10® 

-14.5200 

5 

2.5701 X 10*^ 

2.1300 X lO'* 

3.0110 X 10® 

2.0400 X 10“ 

-3.1691 X 10® 

+0.58265 
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[19] using correlations for the viscosity of the reaction mass and the acti^'-nty coefficients of water 
and caprolactam (see appendix 1). These correlations have been developed [19] by curve-fitting 
one set of industrial data, and have been found to be satisfactory for other operating conditions. 

The simulation package can be used to predict several important properties of the prod- 
uct {Uni PDI, monomer concentration, cyclic dimer concentration, etc.) as well as the various 
reactor characteristics [heat and mass transfer coefficient, vapor discharge rate, V 7 ’(t), temper- 
ature, etc.]. The model of the industrial semi-batch nylon 6 reactor can be used to optimize 
its operation. In this study, we obtain the optimal vapor- release-rate history, V 7 -(t), which 
minimizes the total reaction time, ty, while requiring a few constraints to be met at t=ty. The 
objective function , I, to be minimized is selected as the dimensionless total reaction time 

min /[x(f),t/] = (4) 

where ty^e/ is the total reaction time being used in the reference case (current value). We 
choose the following constraint at the end point: 


?/>! = convtj — convd = 0 (5) 

where the monomer conversion, conv, at the end of the reaction is forced to be equal to a desired 
value, conv^, selected as the current value. This ensures no additional load on downstream 
extraction units. In addition, we use the following stopping condition: 

p = Pn,tf - = 0 ( 6 ) 


where the number average chain length of the product, pri,tp is forced to be equal to a desired 
(current) value, pn,d- In Eqs.4-6, we have (see nomenclature) 


conv = 1 — 


FjC,] + Cl 

Fo[C^]o 


( 7 ) 


pn Pi/ Po 


( 8 ) 
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The stopping condition on ensures the required quality and physical properties of the product 
(the stopping condition and the end point constraints are interchangeable — in fact, if there are 
several requirements at t=t f, any one can be selected as the stopping condition). 

The problem defined in Eqs.4-6 involves three important characteristics of nylon 6 reac- 
tors : monomer conversion. and reaction time. One would also like to minimize the cychc 
dimer concentration, [C 2 ]j^ in the product, simulataneously. This, however, would lead to a 
dynamic multiobjective optimization problem, involving the generation of a Pareto set [16,23]. 
Such a problem is associated with severe computational problems for the semibatch reactor 
studied here, and one must have a good knowledge of solutions of simpler optimization prob- 
lems before one can hope to solve for the Pareto sets. It is hoped that the simpler problem 
defined in Eqs.4-6 will lead to sufficiently low values of [C 2 ]t^ (compared to values encountered 
currently) and will be satisfactory as well as useful for industrial implementation. 

The algorithm for the general optimization problem (with one stopping and m end-point 
constraints as well as n state variable equations and I control variables) [14] is given in appendix 
2. In the present case, three sets of adjoint variables, Af, Xp and Af(f = 1,2, ■ • • ,n) have to be 
evaluated. The ODEs defining these variables are similar in form but their values differ because 
of the different boundary conditions at t=t/. The numerical procedure used is as follows [14]: 

1. An initial control variable history (a single control variable, u = Vr, is used in this study), 
Vj(t), is selected. In order to reduce the memory storage requirements, the values of Vy, 
and the A’s are stored at intervals of 0.1 hr and linear interpolation used. The 15 state 
variable equations with known initial values, x(0), are then integrated in the forward 
direction, using the D02EJF subroutine of the NAG library. The forward integration is 
carried out until the stopping condition of Eq. 6 is satisfied within some tolerance limit 
(10~^ in this study). This fixes the value of tj. During the forward integration, the values 
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of 5f/5Vx and 5f/5x [14] are also computed (numerically) and stored at time intervals 
of 0.1 hr. 

2. The adjoint variable equations are then integrated in the backward direction until t=0, 

using the D02EJF subroutine, and the boundary conditions at t=ty. The modified ad- 
joint variables, and are computed at each storage location (At=0.1 hr). 

In addition, the three integrals, 1$,/ and I//, are computed using the DOIGAF 

subroutine of the NAG library. 

3. Values of two computational parameters, r and e* [14], are selected, and the correction, 

dV'f, computed at each storage location. The new vapor-release-rate history, is 

generated using; 

V^'^^{t) = V-^{t) + 6VT{t) (9) 

4. This completes one iteration. These steps may be repeated (using the new Vx history as 
Vp) until there is no change in the pressure (or Yp) history and the end-point constraints 
are satisfied to a prescribed degree of accuracy (because first order techniques converge 
very slowly as the optimum is reached). 

While carrying out the optimization, it was desired to constrain some of the state variables 
to lie between lower and upper bounds; 

0<n<1.0 (10) 

220‘’C <T< 270°C (11) 

The upper bound on pressure does not allow the pressure rating of the industrial reactor to be 
exceeded while the lower bound ensures the pressure inside to be above atmospheric. The lower 
limit on T is the melting point of nylon 6 while the upper limit represents the approximate 
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boiling point of pure caprolactam (at atmospheric pressures). A standard procedure of satisfy- 
ing these constraints on the state variables, is to incorporate them into the objective function 
as penalty functions. Our previous experience (initial runs) haa shown that the temperature in 
the reactor almost never violates the constraints. Therefore, the only state variable for which 
constraints need to be taken care of, is the pressure inside the reactor. This is controlled by the 
vapor discharge rate, Vx, from the reactor. This suggests an empirical (and easier) approach 
to satisfy the constraints on p(t) in Eq.lO. Some conditional statements have been included in 
the computer code, which take care of the pressure going below or above the lower and upper 
vahies, respectively. Whenever the pressure drops below po, Vx is decreased in infinitesimally 
small steps till the lower bound is satisfied. In a similar way, the upper bound in Eq.lO is taken 
care of. The computer code so developed was found to work satisfactorily. The correctness 
of the code was checked by running it for some special cases. The program was run in the 
simulation mode, using Vx(t) from Ref. 19. The results were in complete agreement with the 
earlier simulation results, thereby suggesting that at least in the simulation subroutine, there 
were no errors. Hand calculations were also performed using the Euler method [22] and some 
of the adjoint variables computed. These matched values generated by the code. These checks 
suggest that the code is relatively free of errors. The optimization code took a CPU time of 
about Ihr 51min on a DEC 3000 ccxp machine (for about 450 iterations). 
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CHAPTER 3 


RESULTS AND DISCUSSIONS 


The reference (current operating conditions indicated by ‘ref’) results axe first generated 
using the simulation code [19] for the three initial water concentrations, [W]o, of 3.45%, 2.52% 
and 4.43% (by weight). The values of several variables describing the current operating con- 
ditions are given in Refs. 19 and 23. These provide values of and in Eqs.5 and 6, 

and also provide results against which optimal solutions can be compared. The current (ref) 
pressure histories for the three values of [W]o are shown in Fig. 2. 

Optimal %/’(!) histories are now generated (with Tj = Tj,re/), using the following param- 
eters in the algorithm (see appendix 2) 


r = 0.1 

(12) 

e* = 0.9 

(13) 


The starting control variable history, V°t, is shown in Fig.3 (case 2). This is selected so as 
to give the starting pressure history shown as Iteration no. 1 in Fig.4. The pressure builds 
up to a maximum value, pmax (< Pmax.ref) in the beginning, is then maintained at this value 
for some period of time, and then decreases slowly to the final value at t=t/ (determined 
by the stopping condition in Eq. 6). Fig.4 shows the variation of the pressure history with 
iteration number during the course of optimization for the [W^jo = 3.45% case. The objective 
function, end point constraint violation, and the pressure history converge (see Figs.4-6) to near 
optimal values in about 332 iterations. In the initial phase of optimization, the convergence is 
observed to be very rapid but then it slows down and some oscillatory behaviour is observed 
before final convergence. This kind of oscillatory behaviour has been observed in earlier studies 
[12-14] also. The convergence of the algorithm can be expedited by manipulating the two 
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f,ref 



Iteration no. 

Fig.5 Convergence of objective function, I (= C/t/,re/) with iteration no. for _ 
case 2 (Fig. 4). 
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Aiion 



Iteration no. 

Fig. 6 i^ii= convt, - convd) vs. iteration no. for case 2 (Fig.4) 


computational parameters, £" and r. In the earlier studies [13,14] it has been established that 
higher values of e tend to speed up convergence. But all our attempts to do so were thwarted 
by the presence of discontinuities in the model. A large value of e* results in large values of 
which lead to large changes in the control vector, VT(t). This results in a sudden drop of 
piessure. The pressure could go down to such an extent that sudden changes in the interfacial 
concentrations of monomer and water are observed. The gradient of the monomer concentration 
(at the bubble surface) then becomes negative. This creates problems in the integration of the 
stat.e- variable equations, since mass transfer from the bubble to the liquid is not allowed for in 
tlie model (and probably does not occur). To avoid such numerical problems, the value of the 
increment, (')V(t), is modulated by a factor u;(< 1). In this study the value of uj has been used 
as 0.75. Therefore, hhpO is replaced by 

I/"'=-(i) = V^{t) + (14) 

in the algorithm. 

One of the most crucial parameters in the control vector optimization technique [17] is the 
choice of the initial control vector history, Vy(t). An improper choice of the initial guess can 
lead to erroneous results or lack of convergence. A common practice followed in most studies 
lias been t,o assume a constant value of the control variable at all values of t. In our study, it 
was not possible to start with such a history since it led to negative gradients of the monomer 
concentration at the bubble surface. Therefore, a Vt° history was selected which corresponded 
to an initial pn-ssure history, p°(t), which was somewhat similar to the one being used currently 
in tlie industrial reactor (see Fig. 2). Optimization was carried out with three different initial 
Vr" liistories (cases 1-3). Fig.7 shows the initial (Iteration No. 1) pressure histories for these 
three cases. The optimal history is found to be almost independent of the starting history 
(see Fig. 8), and the values of t//0.re/ (=0.45) match to two decimal places. Optimization 


27 



Table 2 

Effect Of r and e* On Convergence"^ 
[W], = 3.45% (Case 2) 


e* 

r 

Iterations 

CPU Time 
(sec) 

0.75 

0.50 

374 

6734 

0.90* 

0.10 

332 

6182 

0.50 

0.05 

455 

7088 

0.50 

0.10 

455 

6628 

0.50 

0.25 

455 

7015 

0.50 

0.50 

455 

6946 

0.50 

1.00 

455 

7023 


+ if ft f, ref = ^Ah]convtJconvd = 1.{)Q for all cases 
* bestvalue 

Table 3 

Optimal Conditions for Three Feed Water Concentrations 


% 

If /if, ref 

convtj / convd 

[C2]tjl[C2]tf,ref 

2.52 

0.43125 

0.99956 

0.2526559 

3.45 

0.45000 

1.00000 

0.2584115 

4.43 

0.46250 

0.99977 

0.2491274 
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has also been earned out with different combinations of e* and r. Table 2 gives the results. 

It ca.u 1)0 inferrcKl tliat, tfie final results are independent of the values of the computational 
parajiu'tcu s (pi o\ ideal that the choices are such that negative monomer concentration gradients 
at the l)ubl)l<^ suxlace aie not encountered). It is found that increasing the value of e* leads 
to faster coiiveigeiice (but for certain combinations of the computational parameters, a higher 
valiK* of ('* huids t.o erroneous results and oscillatory behaviour). Although the value of r does 
not affVah, (.he nuniluu' of iterations too much, it has some effect on the CPU time and the 
iuiiDunl. oi oscilliUtjry behaviour. Higher oscillations result in a sudden rise in the objective 
runction which implies more luuiiber of calculations in the forward and backward integration 
,st('ps, thereby result iiii.'; iu higln'r (IPU times. Thus, some amount of numerical experimentation 
lias to be [lerfornied b('furi' t.lu' best choices of r and e* can be deduced. We recommend use of 
values of 1).!) and 0.1 for <* ami r (Eq. 12 and 13), respectively, for the present study. 

Fig. 2 shows t ine optiiual [iressure histories (curves marked ‘opt’) for the three values of 
[W]„. 'rh<‘ [dot. for [W]„=3.4r)% is the same as in Figs.4 and 8. The current and optimal 
vaiior r<dea.s<- raf,<‘K arc shown in Figs. 9-11. Two different maxima in the current Vr histories 
are observed. 'I'he first nia.ximum corresponds to the opening of the valve for the first time 
to a.rr<'st, th(> increasi' iu iiressure associated with early vaporization, while the later one is 
aHsoeiahed with the desired ra,[)id fall of the pressure. It is evident from Figs.9-11 that different 
Iiressure liistori.'s ami vapor release rates are required to produce three different grades of 
{lolyiner. Table 3 give.s t,he oiitimal conditions for the three water concentrations. 

Fig, 12 shows the variation of /q. with t for the optimal as well as the reference runs for 
(,lu- three [W]., values. In tlu« r(d-erence (current) runs, two distinct regions are observed where 
tlie lirst, being associated with the polyaddition reaction while the second with 
tl... n.,,d,iou. Tl,erc is a short plateau m-between the two regimes. From the 
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vj ' vj max, ref 



Fig. 9 Dimensionless vapor release rate histories. Curve (ref): reference run, (opt) 

optimal run, for Tj = Tj,re/, [W]a = 3.45%. 



VT/Vj^max/ef 



Fig. 10 Dimensionless vapor release rate histories. Curve (ref): reference run, (opt): 
optimal run, for Tj = Tj.re/, [W]o = 2.52%. 
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^ ^T, max, ref 



Fig. 11 Dimensionless vapor release rate histories. Curve (ref); reference run, (opt): 
optimal run, for Tj = [W]o = 4.43%. 
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Fig. 12 Variation of degree of polymerization with dimensionless time for the 
reference (ref) and optimal (opt) runs for Tj = Tj,rej, [W]o = 2.52%, 
3.45%. 4.43%. 
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figure it is evident that this plateau does not exist for the optimal histories because of early and 
rapid pressure decrease, and the increases rapidly at almost a constant rate to the desired 
value. The value of approaches its final value almost asymptotically at t j^ref for the reference 
runs. This is not so for the optimal runs. This makes it necessary to have an excellent control 
for the reactor so that the value of fj.^ of the product does not exceed the desired value during 
emptying of the reactor at the end of the run. With the availability of highly sophisticated and 
robust controllers, this should not pose any major difficulty. 

Fig. 13 shows the dimensionless temperature histories for the reference runs as well as for 
the optimal conditions. In the initial phase, when pressure in the reactor is building up, the 
temperature variations for the reference runs are similar to those for the optimal runs. This is 
not surprising since the pressure histories for the optimal and reference runs are quite similar 
(see Fig. 2) in the beginning. However, the temperature histories for the optimal cases deviate 
from the reference values when the optimal pressure histories start deviating from currently 
encountered values. This suggests the importance of latent heat effects. Fig. 14 shows the 
variation of monomer conversion with t, for [W]o=3.45%. It is observed that the conversion 
rises less sharply for the optimal conditions, than in the reference case. This is attributed to 
lower temperatures in the optimal case. 

Fig. 15 shows the variation of the cyclic dimer concentration, [Ca], with time. The cyclic 
dimer concentration remains negligible upto some time and then it starts increasing. The 
final value for the optimal case are found to be considerably lower than the current values. 
This is a blessing, since we had obtained the optimal Vr(t) without putting any requirements 
on the cyclic dimer concentrations. It must be kept in mind that the low concentration of 
[CJ,, for optimal runs could increase further during post-polymerization processes (including 
the emptying of the reactor), since equilibrium is not attained for the reactions involving the 





PaUOO / AUOQ 



time for the reference (ref) and optimal (opt) runs for Tj = 
[W], = 3.45% (case 2). 



[C2] /IC 2 I tf>ref 



^ ^ ^f/ef 


Fig. 15 Variation of the dimensionless cyclic dimer concentration with dimensionless 
time for the reference (ref) and optimal (opt) runs for Tj = 

[W]„ = 2.52%, 3.45%, 4.43%. 
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cyclic dimer. However, since one would be ensuring that the values of /i„ do not go up beyond 
the desired values during emptying of the reactor, this requirement should pose no additional 
problems. 

It is interesting to note that the optimal pressure histories obtained (at Tj = Tyre/) in 
the present study using Pontryagin's minimum principle, is quite similar to that obtained by 
our group earlier [23,24], in which the shape of the pressure history was fixed a priori, and 
was determined completely by a few constants {prnax^ tc; S, ty). Optimal values of three of 
these constants (p^naer: 0 and S. with tj determined by the stopping condition of Eq. 6) were 
determiined using sequential quadratic programming. The similarity of the optimal pressure 
histories as obtained more rigorously in the present study, with the optimal histories obtained 
with somewhat artificial constraints put on the shape of p(t), confirms the validity^ of the 
constraints used in our earlier study. It may be added that major increases in the plant 
capacity have been achieved, using some of the results of the present study. 

Fig. 16 shows how the optimal pressure history changes when the jacket fluid temperature 
is increased by o'^C. Higher pressures are required (to suppress the vaporization of monomer 
and water). Use of higher values of Tj lead to lower values of the objective function. However, 
it must be ensured that the polymer does not get degraded under these conditions. 

In addition to the optimization problem studied above, we also studied the following 
problem (for [W]o=3.45%) 


min /[Ur(0,Tj = 


with the following end-point constraint/stopping condition 


(15) 


[C2],J[C2]t,,reJ = 0.2584115 


(16) 







1-20 


[Wlo =3-45y. 


1-0 


0-80 



Fig. 17 Optimal pressure histories for the two optimization problems described 


in Eqs. 4-6 and Eqs. 15-17 for [W]o = 3.45 %. 
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The values of and conv,, /conva under optimal conditions came out as 0.4375 and 

0.96689 respectively, which compare very well with the corresponding values of 0.43 and I.OOOO 
for the problem described in Eqs. 4-6. The opt.mal pressure histories are shown for the two 
problems in Fig.17. The results for the two problems seem to be in close agreement and indicate 
that the optimal results are not affected, significantly, by the type of end-point constraints used. 



CHAPTER 4 


CONCLUSIONS 

Optimal vapor release rate histories for an industrial semi-batch nylon 6 reactor, using 
Pontryagin's minimum principle have been obtained. The results predict a significant reduction 
in the final reaction time while sim.uitaneously (and somewhat fortuitously) leading to fairly 
low values of the cyclic dimer concentration in the product. These optimal histories can be 
implemented on the existing reactor without major alterations and investments. The compu- 
tational skills developed in this work can also be helpful in the solution of more complicated 
optimization problems, e.g.. multiobjective optimization. 


44 



CHAPTER 5 


SUGGESTIONS FOR FUTURE WORK 


Though this work has been able to generate optimal pressure histories which could lead 

'h 

to substantial amount of improvement in the performance of the industrial semi-batch nylon 6 
reactor ~ interms of total reaction time as well as concentration of undesired cyclic dimer, still 
it needs to go some steps before it could attain the culmination. 

In this study we have assumed that the other manipulated variable i.e. jacket temperature, 
Tj, is kept constant throughout the reaction. A more realistic solution of the problem would 
have been the one in which this variable is included as second control variable. The two control 
variable problem, which is computationally intensive, was also tried. The convergence was 
found to be very fast, but the discontinuties in the model did not allow the optimization to 
proceed after some itcirations. Therefore it becomes neccessary to remove these discontinuities. 

The inclusion of second end-point constraint, dimer concentration [C 2 ]tj^ will lead to 
multiobjective optimization problem. The solution of this problem is called for so that pareto 
sets are generated to chose the reactor operating conditions as per the need. 
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Appendix 1 


Mass and energy balance equations 

^=^.CC.][W]-J:;cs.]-2fc,[s,K+2;:;[\vi(fi,-[S.])-A-s[S.][C,]+^-s[S2]-/c5[s.][c,]+fc5[S3]-[s,]— — 

^ = fc, [C J[W] - fc; [S j - ic, uf + fc;[\v](;,, + fc^[\vi[C,] -A-;[S,] +/i„ 

^ = t:.[C,](\V]-^:;[S,] + <:3[C.3;t„-A.--3(M.-[Si3) + 2A:,[C,];.„-2r.Ut.-[S,]-[S,]) + 2fc^[W][C,3-2A;[SJ-rM, 

iiu-. 1 

-^ = A-i[CJ[W3-fc3[S,3 + 2i3pi-r-*;[W'3(;i,-P3) + fc3[C.3(;t„ + 2,t,) + /c3(;i„-2/t,+[S.3) + 4A'5[C33(;i,-^,|+4A-3(,i3-;i, + [S33l 


+4fc.[W][C33-4A;;[S33+M3 


0.n3J?„„-0.018R,. 


^ = - t.[C33 W + k'.lS,-] - fc3[C33/^o - - [S .] - [S 33 ) + [C 33 

^=-A-.[C.3m+t3[S.3+fc3;^o--^-.[\V3(u,-/ij-ii.-.[C33[W]+^’[S,3-Aw/f+[W3 °’ - ’^^~"'^°'°‘^^'’^' 

df 

— =-{0.I!3K.„, + 0.018f:,J 

dr 

dT f F 5 ') 

_=|{..U(rj--n + — £ r,.t-AH,)-[R,„;.„(r,) + R..;.jr,)3-[0.113R.„C;,„4-0.018R..«C,.J(r-T,)^CU.C0.113K.„-r0.018R.J(T- 
>= {[CUx + 10925 X 10- ^{T- r,)3R}- ‘ 



iV[M'3 


dt ■ 1'. 

K,([M'3 + [W''3 + 

[N'll 

d[W'] R.. 

i't[\V''3 


dt ■ 1'. 

K,([M''3 + [W'3 + 

[N']) 

d[N«]_ 

P'tCN'I 


dt i; 

,([\r3 + [w'3 + [N' 

1) 


Closure condiiions: 
[S3] = [S,] = [SJ 

^3 = 
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Correlations and other equations used 


Expression for Vj 
Stage 1: 

Stages 2-5: = 


RT\dt 


[M'']4-[W''] + [N'] fdT\ 
' T V"^/ 


Rates of vaporization 

-[C,],*) [C,]r = [M-'isr 

i?v. = f(fc,..a)f([W] - [W]*) + - Ol;) 


Pressure 

P=;[Nr]-f.[WT + [N^];jRr 


Equations for activity coefficients 

f[C,] 


monomer conversion = 1.0 — 


foCc.i-:. 

,, . /^mf ““ l^ma , 

, = /Lo-l (monomer conversion) 

0.95 

i ^ w i w o 

. = H (monomer conversion) 

0.95 


1 apnur pressure 

7024.023 

ln[P;"‘(kPa) 101.3] = 13.0063 (rci. 16i 

3816.44 

InfP'J' (kPa), 101.3] = 11.6703 (ret. IT) 

T(K.)-46.13 


Diffusion coefficients 
(m- h’M = 3.6x 10"'’ 

£?„{m“h-') = 2.88xlO-*' 

Latent heats of vaporization 
r, = 473.15 K 

r,) = 34.2559 kJ mol"' (ref. 18) 
;.jr,) = 51.0193kJmol" ' (ref. 18) 
C;.„ = 1.6426 kJ kg-' K"' (ref. 18) 
C;., = 1.9963 kJ kg-' K"' (ref. 18) 


Heat transfer coefficient 


fq (kJ h"' m“- K-') = 


[//(Pas)]'^ 


(7=- 


l thickness (m) 
-H 


Correlations from ref 19 for mixture physical properties (liquid) 
(kJ kg"' K'‘) = 2.0925 + 2.0925 x 10"^ [T(1C)-273.15] 
p (kg m"'')= la'lO;!. 1238-0.5663 X 10"^ [T(K)-273.15]J 
k (kJ h"' m"' K-') = 0.7558 


Correlations from ref. 20 
d^np 






6;7( poise) 
360//( poise) 


P^i 

[C,] 


[C.]+[W] 


.. p«. 

; m-‘ m 


for mass transfer coefficients 


Lcj-tyo 



(i) Quiescent (q) desorption iQ^<P) 

i = m or w 






_ ikl.i)(.ci^r 

£?, 

{nj4)D; 

(f/p) 


(^■l.w^)h==0 


(ii) Bubbly (b) desorption ifL, > 
_ [W]-[WJ^ 

[Wh]" 

cr, = I.81iVR-,‘^*^ 


If < <7,.: 

{^Lw^)h {h’‘) = 6.77 X 10■^^•;;,-T-•"'' X 3600 
1^1. If 


If fT><7c: 

(ic,.^a)b(h-‘) = 2.45 x x 3600 

log,o0 = 522((j-(7jNR-;«' 
ik\,m^)f ~iky ^a^f ^^(j) 


Kiscastty correlation {neglecting effect of wafer) 


Ucp) = 2.7969 X 10"^ exp[3636.364-T{K)] for r>473.15 K 
V kg polymer J c-o «■' c-o rjm c 


1875.0 

riK) 


^4.678 



0.75 


M„=113 Mw=H3 

c (kg polymer/ 100 kg mixture) = 1 1.3 


For c<2i.0/M 

-!?^=1.0-0.3102cM+0.0575(cM)"-0.525xl0--(cM)^+0.r-05 

X 10'^(cMr -0.3663 X 10-’(cM)= 


For O21.0/M 

log,o >1 (poisei = 5 log, 12.3097 
for C;.\/““*>315 and .\/.>5000 
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else 

logio q (poise) = log^olQ^f v.) - 3.503 

where C' (g polymer/cm^ mixture) = 11.3 x i0~^p/ij 



Appendix 2 

EQUATIONS FOR OPTIMALITY WITH END POINT CONSTRAINTS AND 

STOPPING CONDITIONS 


1. Objective Function 


min 

2. Mass Balance (State Variable) Equations 



■ dx] ‘ 
dt 




dx 

dt 

dxr) 

dt 

= f[ 

)rU(t) ] = 

f2{x,u) 


dXn 

- dt - 



_ fn{x,u) _ 


with 


u^(t) = [ ui, U2. • • • . U/ ] 

3. End Point Constraints 


^[x(^/)V/] 


Vi(x(t/)V/) 

Vm(x(t/)V/) 


4. Stopping Conditions 


d[x(^/), tf] = P • • • , tf] = 0 

5. Constraints on Control Variable u 


Umin < U < Ur 


6. Adjoint Function Equations 

(a) 


d\^(t) d 

dt dt 


A' 


Xi 


n J 




(—] 


where 
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1.1. f.. lUWIfW 


A^(t) 




■ dfi df2 

.iln 1 

dxi dxt 

dxi 

^ili . . 

. ^ 

8x2 8x2 

8x2 

dfi 8/2 

dfn 

. dXn dXn 

dXn . 



(b) 


dX^jt) _ d 
dt dt 


A Vl \ ^2 

^2 A2 


Aj’'" 






-A^(t)A^(i) ; X^(tf) = 


'd^\ 

. j. 




:c 


dX^(t) d 

dt dt 


Increment In u(t) 


Af 

^2 




^ -A^{t)XP(Vj : X^(tf) 


.dy:)t=tf 


du(t)^ui) = 

+ 


X 


where 




I//(lxl) ^'®'5(mxm)^^I(”^x 1) 

Bflxn)WAf„xm)Wl$^^(mxm)^^ 


(mx 1) 


T P 

^^^{mxm) — Jq y^{mxTi) 




t/;,lxl) = Id (^{Ixn) 


= ^(nxl)(*) 


T 


B(„,i)(i)Bfi,„)(i)A®/,„(t)dt 

1 

B(nxl)(^) Bflxn)WAjLl)(t)ht 

B(nxi)WBfi,<„)(i)Ajj^i)(t)dt 


A(nxl)('0 


Qi ^ \dx) ^ 

§£. _L f^V'^f 

5< ^ VaxJ ^H=i, 
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B 


T _ 

Ixn 



MlUz. . 

. . M 

dui dui 

dui 


. . M 

dun du2 

du2 

MM. 

..iln 

du] dui 

dui 


8 . Increment 


(5^(^xi = -c*^(t/);0<c<l 


T ’ 

Ixm 

f 

t=t j 



